Thermal conductivity for a chain of anharmonic 
oscillators perturbed by a conservative noise 



Cedric Bernardin 



Abstract We consider tZ-dimensional chains of (an)harmonic oscillators we perturb 
by a noise conserving energy or energy and momentum. We review the thermal 
conduction properties we obtained for these systems (|[8l. HI, 01, 0) and conclude 
by several open questions. 



1 Introduction 

The derivation of phenomenological laws from a microscopic description of the 
matter is one of the goals of statistical mechanics. Among them Fourier's law is 
probably one of the simple: when a small gradient VT of temperature is applied to 
a material, in the steady state, the energy flux J is proportional to the gradient of 
temperature 

J = -kVT 

The proportionality coefficient K is called the thermal conductivity. Despite its sim- 
plicity and the interest it has in the physical and mathematical community the deriva- 
tion of Fourier's law from a microscopic model remains one of the main open ques- 
tion of nonequilibrium statistical mechanics (" 111711 . ifTTl . Q3]). 

In insulating crystals heat is transported by lattice vibrations, and since the pio- 
neering work of Debye, systems of coupled anharmonic oscillators have been used 
as microscopic models for heat conduction. They are classical system of particles 
interacting through a nearest neighbour interaction potential V and which are in an 
external potential W. The Hamiltonian Jt? is given by 
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* = I ( fep + + £ v fe - «*) J > A c z " 

aga \ x b-*l=i / 

where m x ,q x ,p x denotes the masse, position and momentum of the particle with 
equilibrium position x E A. 

It is well known that harmonic chains, because of their infinitely many conserved 
quantities, have infinite conductivity and do not obey Fourier's law ((T9|)- This be- 
cause phonons can traverse ballistically along the chain. It is often expected that 
enough strong nonlinearity or disorder (like the presence of random masses) causes 
scattering between phonons and should imply a sufficiently fast decay of correla- 
tions for heat current and hence a normal conductivity. A rigorous treatment of a 
nonlinear system, even the proof of the existence of the conductivity coefficient, 
seems to be out of reach of current mathematical techniques. In this context the 
understanding of the coupled effect of nonlinearity and disorder is a challenge. 

The situation is in fact more complex. In some low dimensional systems (d < 2) 
anomalous thermal conductivity is observed numerically and experimentally in nan- 
otubes technology. The anomalous conductivity in low dimension has attracted a lot 
of attention in the literature, and it has been suggested that conservation of momen- 
tum is an important ingredient. There is no agreement, theoretically and numerically, 
about the exact dependance of the conductivity with the size of the system ( ifTTl ). 

Hence it makes sense to look at simple models which incorporate the important 
features that one believes are necessary to see normal transport. The main difficulty 
in Hamiltonian dynamics with a large number of degrees of freedom is to show 
that they behave ergodically, e.g. that the only time invariant measures locally ab- 
solutely continuous w.r.t. Lebesgue measure are, for infinitely extended spatial uni- 
form systems, of the Gibbs type. For some stochastic lattice gases it can be proven 
but it remains a challenging problem for Hamiltonian dynamics. Taking advantage 
of mathematical techniques developed in the hydrodynamic limits communauty we 
introduce hybrid models between purely Hamiltonian systems and purely stochastic 
models which remain mathematically tractable but are sufficiently close to realis- 
tic systems to reproduce at least qualitatively what is observed for these systems. 
We consider chains of oscillators perturbed by a stochastic noise conserving en- 
ergy or energy and momentum. These stochastic perturbations are here to simulate 
(qualitatively) the effective (deterministic) ergodicity coming from the Hamiltonian 
dynamics. 

The paper is organized as follows. In section|2]we introduce a model of coupled 
oscillators perturbed by a noise: the first noise conserves only energy and the second 
one conserves energy and momentum. In section[3]we review linear response theory 
and Green-Kubo formula. Section|4]is devoted to the study of a chain of oscillators 
with the noise conserving only energy. For the harmonic homogenous harmonic 
chain we show Fourier's law is valid and compute the conductivity. For anharmonic 
chains we provide lower and upper bounds. The effect of disorder is considered 
in subsection 14.31 Then we consider the energy-momentum conserving model in 
section|5]and show that in the homogenous harmonic case a breakdown of Fourier's 
law holds for low dimensional momentum conserving systems. We provide also 
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upper bounds for the conductivity in the anharmonic case. We conclude the paper 
by open question in section [6] 

Notations : The canonical basis of W 1 is noted (ci,C2j • • • , e d) and the coordinates 
of a vector u G W 1 are noted (a 1 ,. . . ,u d ). Its Euclidian norm \u\ is given by \u\ = 

(u 1 ) 2 + . . . + (u d ) 2 and the scalar product of u and v is u ■ v. 

If N is a positive integer, denotes the cZ-dimensional discrete torus of length 
N and we identify x = x + kNej for any j= l,...,d and ieZ. 

If F is a function from If 1 (or T^) into R then the (discrete) gradient of F in the 
direction e ; - is defined by (W ej F)(x) = F(x + ej) — F(x) and the Laplacian of F is 

given by (AF) (x) = Z d j=1 {F{x + ej ) +F(x -ej) - 2F(x)}. 



2 The models 

In this section we introduce deterministic nonlinear chains on a multidimensional 
lattice perturbed by a stochastic noise. The stochastic perturbations are such that 
they exchange momentum between particles with a local random mechanism that 
conserves total energy or total energy and total momentum. 



2.1 Closed system 

We first consider the closed system with periodic boundary conditions. The atoms 
are labeled by x <E Tjy. Momentum of atom x is p x € R , its displacement from its 
equilibrium position is q x e W 1 and its mass is m x > 0. The configuration space is 
given by Q N = {(q x ,p x ) € R (/ X W 1 ; x G T^}. The Hamiltonian is given by 

We assume that V and W have the following form: 

d d 
V(q x ~ qy) = E Vj(qi - q{), W(q x ) = £ Wy(^). 

and that V;,Wi are smooth, non-negative and even. We call V the interaction poten- 
tial, and W the pinning potential. The case where W = will be called unpinned. 
In the sequel we will refer to the (a, v)-harmonic case: 



2m x 



■W(q x )- 



E V{q x -q y ) 

y-x\=l 



Vj(r) = ar 2 , W j (q)^Vq 2 1 a > 0, V>0 



(1) 
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for which explicit and generic results can be obtained. 
The generator Jzf/v of the dynamics is defined by 

-SSv = S*N + Y^N 

where s^n is the Liouville operator corresponding to the Hamiltonian 3% 

= E { d Px-#& ■ d gx - d qx J4? N ■ d px } 

X 

and is the generator of the Markovian noise. The parameter 7 > regulates 
the strength of the noise. It acts only on momenta and is local. It consists in an 
infinitesimal exchange of momenta preserving some conservation laws. The first 
conservation law we impose is the energy conservation. Since the noise acts only on 
momenta it is equivalent to require conservation of kinetic energy. This corresponds 
to the so-called energy conserving noise. If we require also the conservation of total 
momentum then we get a second noise we call the energy-momentum conserving 
noise. Let us denote % x = p x j \fm~ x . In the homogeneous case m x = 1, % x = p x . 
The two noises have the following form 

**=\t E ra 2 

\x-z\=l 

where X' x l is equal to 

for the energy conserving noise and equal to 

K: = i4- - d 4 ) - (nt-nbfy - d^. 

for the energy-momentum noise if d > 2. If d = 1 in order to conserve total momen- 
tum and total kinetic energy, we have to consider a random exchange of momentum 
between three consecutive atoms ([2]). 

The interpretation of the vector fields is the following. To be specific we take d = 
1 and consider the energy conserving noise. In this case X x =X C ' j = (n x+ id ny — 
% x d% x ). Observe thatX x is the vector field tangent to the circle C x = {(n x , TTx+i); 7t x + 
K x+l = 1} so that generates a diffusion on C x . In fact it is nothing else than a 
standard Brownian motion on the circle. The generator S^n corresponds to a sys- 
tem of coupled Brownian diffusions preserving the kinetic energy L vG t<' W 2 - The 
energy-momentum conserving noise is defined by a similar procedure but the sur- 
face C x has to be replaced by the surface of constant kinetic energy and constant 
momentum. In dimension 1, this surface is reduced to a point and it explains why 
we have to consider a three-body interaction. 

Because the noise conserves energy, a family of stationary translations invari- 
ant probability measures for S£n is given by the Gibbs measures. In the energy 



Thermal conductivity for a perturbed nonlinear chain 5 

conserving case they are parametrized by inverse temperature /3 = 7" 1 and in the 
energy-momentum conserving model by inverse temperature j3 and mean momen- 
tum average p. We denote the Gibbs measure with inverse temperature j3 = T~ l and 
zero momentum average by ji N T . It is given by 

H N , T (dqdp) = Z^exp{-^J%)dqdp 

where Z N T is the partition function. Expectation with respect to fi N T is denoted by 
(■)n,t- Remark that in h 2 (jj, NT ) the Hamiltonian vector field stf^ is antisymmetric 
and the noise is symmetric. 



2.2 Open system 

We now consider the case where the system is in contact with thermal baths at 
different temperatures Tf and T r . Thermal baths are given by Ornstein-Uhlenbeck 
processes with the corresponding temperature. To simplify notations we take d = 1. 
The configuration space is now given by Xn — { (Px,1x) 6lxl;i=l,... ,N}. The 
generator of the evolution has the form 

N 

£n = £ {d Pz J0h ■ d qx - d lh Mk ■ d Px } 
.v=i 




We have to specify boundary conditions for MU- For example one can define 

iV-l 
x=l 

with qo and qN+i fixed. 

Even if one can prove there exists a unique stationary probability measure 
for the process, there is in general no formula to express it. The only case where 
one knows (•)«.„ is the equilibrium case 7) =T r = T where the stationary measure 
is given by the Gibbs measure at temperature T. Otherwise the probability measure 
(•)jv,.„ is called a nonequilibrium stationary state. 



2m x 



W(q x ) + - £ V(q x -q y ) 



3 Thermal conductivity and linear response theory 



In this section we review briefly linear response theory and Green-Kubo formula for 
thermal conductivity. Derivation of the Green-Kubo formula is heuristic and even 
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its (mathematical) existence is a challenging problem. Roughly speaking Green- 
Kubo formula is the space-time variance of the total current at equilibrium. Thermal 
conductivity is a transport coefficient defined by considering the system out of equi- 
librium. Linear response theory express the fact that if temperatures T( and T r are 
different but close then a linear approximation is valid and the thermal conductiv- 
ity k(T) is equal to the Green-Kubo formula K GK (T). Such a formula belongs to 
the family of fluctuation-dissipation theorems since it relates dissipation (i.e. the 
thermal conductivity) to fluctuations (i.e. fluctuations of the total current). 
Defining the energy of the atom x as 

&=^-p 2 x + W(q x ) + - £ V(q y -q x ) 

x r-\y-*\=i 

the energy conservation law can be read locally as 

k=\ 

where J X}X+e ,([Q,t]) is the total energy current between x and x + e^ up to time t. 
This can be written as 



Jx,x+e k 



([0,t])=J j Xl x+e k (*) ds + M x ,x+e k (t ) 



In the above M x>x +e k (0 are martingales that can be written explicitly as Ito stochastic 
integrals. 

The instantaneous energy currents j x .x+e k satisfy the equation 

d 

■^N&x = 52 {jx-e k .x —j Xr x+e k ) 
k=\ 

and it can be written as 

jxa+e k = j x ,x+e k + Yjx,x+e k (2) 

The first term in (f2]) is the Hamiltonian contribution to the energy current 



fL+e k = -\(yV){q x+ e k ~ q x ) ■ ( — 

L \ m x+e k ™>x / 



1 ^ ... , „ / ^ 



7 y J\ ( ix+e k ^x) I y/ 2 ^ 1/2 
7=1 \ m x+e k ™x / 



while f x x+e is the noise contribution. 
For the energy conserving noise we have 
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1 , 

fx,x+e k =-jVeM*\ 

In the energy-momentum conserving case, in d > 2, it is 

fx,x+e k = -Ve k \7tx\ 2 

and in d = 1 is 

Jx,x+i =-V<p{n x -i,nx,%x+i) 
<p(n x -i,n Xi n x+ \) = gfal+i + 47t x+ n x-i + ^+1^-1 — 2^+1^-2^^-1] 

The particular form of j s is not very important. What is relevant is the fact that 
Jx x+et ^ s a discrete gradient V ek ^x of a local function <j> x . It means a weak form of 
Fourier's law is valid at the microscopic level as soon as we can relate (j) x to the 
(macroscopic) local temperature. Observe that it is not the case for the Hamiltonian 
part of the current and one of the main difficulties is to express j" as the sum of a 
discrete gradient and a small term (see ([Toll). 



3.1 Nonequilibrium setting 

We consider the system out of equilibrium (see subsection l2.2b in contact with two 
heat baths at different temperatures T( and T r in the first direction. The conductivity 
k(T) is defined by the thermodynamic limit 

,rr\ r ,. (jO,ei)iV,sj 

KIT) = lim hm 

V ' T(,T r —*T N—>°° Ti>-T r 

We have seen in the introduction one expects that for one and two dimensional 
systems conserving momentum such a limit is equal to infinity. In order to estimate 
this divergence one can study the finite size thermal conductivity 

Kn(T)= hm (3) 

T e ,T r ^T lg~l r 

Under suitable conditions his quantity is well defined even for purely Hamiltonian 
chains but it is not straightforward. It is often expected that 

lim k n (T) = k(T) 



but it is not obvious since there is an exchange of limits. 
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3.2 Green-Kubo formula 



The difficulty arising in the study of k(T) is that we have no explicit representation 
of the nonequilibrium stationary state (-) N . SS . Since T( and T, are close it is suggestive 
to use a perturbative approach to compute k(T). 

Performing a first order development in the stationary state (-) N , SS we get that 

K(T) = K GK {T) 

where the Green-Kubo formula for the conductivity K GK (T) is given by ( ll23l . pp. 
188-190) 



K GK (T) = lim fim ^ £ E N T ([0,t])J 0>ei ([0,/])] (4) 



Here Ea, r indicates the expectation with respect to the equilibrium dynamics starting 
with the Gibbs measure < • > N T at temperature T. This definition itself is formal 
since we have to prove existence of the limits. 

By standard stochastic calculus and a time-reversal argument (O) one can es- 
tablish the following equality 



1 



2T 2 t 



£E^(/^ +ei ([0,f])/ 0iei ([0,f])) 



(2T z N a t) 



i 



(s)ds 



1 

d 



(5) 



Here the term y/d is due to the presence of the noise. This is the first term 
which is of interest. In view of the Green-Kubo formula the anomalous behavior 
of the conductivity should appear in a slow time-decay of the time correlation of the 
(Hamiltonian part) of the current. 

Let us denote 

xeY» N 

In order to study the large time behavior of 



C{t) = lim C N (t), C N (t) = {2T 2 tN d Y l E NT 



we study the asymptotics as — * °° and then A — > of the Laplace transform £at(A) 
of tC N {t) 



£*(A) 



tC N (t)dt 



By stationarity and integration by parts, we have 
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W) = 



1 



Ae-*%, r [3 ei (03ei(0)] 



Denote by e'-^ the semigroup generated by and remark that 3 fl if) = e'^ N Ze x 
then 

1 ^ _j 
£jy(A) = ^ 2 ^ 2 (3 gl , (A — Jf N ) Zei/N,T (6) 

A normal finite conductivity corresponds (in a Tauberian sense) to a positive 
finite limit of A 2 £w(A) as /V — » °° and then A — * 0. In this case, the conductivity 
k(T) should be equal to the following form of the Green-Kubo formula 



K GK {T) = T/d+ Hm lim N~ d (3 ei , (A --2fr) -1 3 ei )w,r 



3,-»0iV- 



(7) 



No general argument gives the existence of the limits in (01 and in (|7]i nor that if 
they exist they are equal. 

Observe that Green-Kubo formula (|4]i predicts only the value of the thermal con- 
ductivity k(T) (defined in the nonequilibrium setting). If the thermal conductivity 
is infinite it says a priori nothing about the behavior of the finite size thermal con- 
ductivity KJv(r) defined by (01. To overcome this problem we define the truncated 
Green Kubo formula by 



Kg K (T) = (2T z N a t N )-% 



\-1ti 



',v 



1 

d 



where t^ = N/v s with v s the sound velocity defined by 



lim 



co(k)\ 



and where 



1/2 



a{k) 



E^/(o)- 

.7=1 



-4y/'(0)sin 2 (7r^') 



is the dispersion relation of the approximated linear system. This definition of the 
conductivity of the finite system is motivated by the following consideration: in 
the harmonic case the finite size thermal conductivity K)v(T') can be obtained by 
this truncation technique (in a rigorous way) and we expect this is still valid for 
the anharmonic chain. In the linear interactions approximation V^ffl^) is the group 
velocity of the £-mode waves, which are the heat carriers, and typically v s is an upper 
bound for these velocities. Consequently tN is the typical time a low k (acoustic) 
mode takes to cross around the system once (see(T7|). Typically v s is of order one 
and we will take v s = 1 in the sequel. 
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4 Energy conserving model 

In this section we state the results obtained for the energy conserving noise. As 
expected one can not obtain infinite conductivity for this model since momentum 
is not conserved. In the homogenous harmonic case one computes explicitly the 
Green-Kubo formula K GK {T) and also the conductivity k(T) in the one-dimensional 
unpinned system. It turns out that k(T) = K GK (T) so that predictions of linear re- 
sponse theory is valid although the proof does non use perturbative arguments. In 
the anharmonic homogenous case one can establish lower and upper bounds for 
K GK (T) indicating a positive finite conductivity. The proof of the convergence of 
the Green-Kubo formula in the anharmonic case remains open. We are also inter- 
ested in the effect of disorder (random masses) in the energy conserving model. This 
is the contain of subsection l4.3l 



4.1 Homogenous harmonic chain 

Here we consider the homogenous (a, v)-harmonic case ([U with all masses equal 
to 1. 

Let us define 

r ( 4a 2 lf_,sin 2 (7r^') \ . . 
D= / ^— \dt x dt d 

We have the following theorem 

Theorem 1. K GK (T) defined by (@ is finite (pinned or unpinned) in any dimension 
and given by 

rr D y 
k gk (T) = -+- j 
y d 

Proof. Recall ©. A simple but crucial computation shows that 

(A— _S?Af) l Ze 



A + 7 



Let Dn = £>A?(a, v) be the constant 

One computes easily limN^oo(^ ei ,Z ei )N.T and after inversion of the Laplace trans- 
form, one gets 

C„(,) = ^(l + I(l-^)) 
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Then theorem follows. 

Observe that if the noise becomes weaker (i.e. 7 — > 0), we obtain a purely ho- 
mogenous harmonic chain and the thermal conductivity is infinite. 

In the following theorem we study the one-dimensional unpinned system in con- 
tact with two heat baths at temperature T( and T r . We show that conductivity is finite 
and coincides with the Green-Kubo formula. The proof remains valid in any di- 
mension and with or without pinning as soon as we are able to prove the following 
bound 

Vxe{l,...,N}, (e x } NiSS <C (8) 

where C > is independent of N. Unfortunately it has been proved only in the one 
dimensional unpinned case. 

Theorem 2. $8$ Consider the one dimensional (a, v)-harmonic case (0 with v = 0. 
For any 7 > 

Urn JV < j x , x+l >= a (7+ 7- 1 ) (2i - T r ) (9) 

Hence we have k(T) = K GK (T) = a(y+ y~ l ) 

The proof of this theorem can be found in |8). It is based on entropy production 
bounds and use an explicit decomposition of the current j xx +i as the sum 

A..v+i =V0 A +JVi, (10) 

where h x , <p x are two explicit local functions and V is the discrete gradient. Hence 
the current is the sum of a dissipative part (a spatial gradient) and a fluctuating part 
(a time derivative). For this reason we call this equation a microscopic fluctuation- 
dissipation relation (0,0>0)- 



4.2 Homogenous anharmonic chain 

The introduction of nonlinearity in the Hamiltonian dynamics complicates consid- 
erably the problem. We do not have any proof of the existence of K GK (T) nor k(T). 
At least we have estimates which indicate that a finite strictly positive conductivity 
is expected. By strictly positive we mean that the Hamiltonian contribution to the 
Green Kubo formula (the second term in (|7])) is strictly positive, the term y/d being 
only due to the noise and of no interest for the study of conduction properties of the 
underlying deterministic dynamics. We consider the (formal) Green-Kubo formula 
©. 

Proposition 1. Assume that 

r''(3 C| ,3 fI ) w <c 

with a constant C independent of N. Then there exists a positive constant C' inde- 
pendent ofN such that for any X > and N, 
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N- d (Z ei ,(X-^ N )-^ ei } NtT <C (11) 

Observe that the assumption done in this proposition is natural and is satisfied 
for reasonable potentials V,W (see (fT2]i). This shows that if Green-Kubo formula 
(0 converges then K GK (T) is finite. 

Proof. Introduce the following so-called H\x and H~\\ norms defined by 

ll/llluH/^-y^/kr 

We have the following variational formula (J6)) 

(&i , (A - j2V) -1 3ei)iV,r = sup <2(Z ei ,u) NJ - \\u\\\ X - \\^NU\\ 2 _ lX } 

To obtain the upper bound we forget the term ||^r«||lj j and we get 

(3ei , (A - -2V) -1 3ei)tf,r < sup | 2(^ ei ,M) A r ir ~ IMIl A } 
It is well known that 

SUp \2(3 ei ,u) NJ - \\u\\ 2 Ll \ = (3et , (A - Y-VNy^Ze^N.T 

and a simple computation shows that 

(A-yj?5v)-% 



A + y 

We are left to prove that 

N- d (Z ei ,Ze,)N,T<C 

with a constant C independent of N. It is exactly our assumption. 
We expect also a lower bound of the form 

C~ x < N- d {Z ei ,{X-X N )- l Z ei )N,T 

with C > independent of X and AT. 

The strategy to prove such a lower bound is straightforward: we have to find a 
good test function Vjv,a and show that 

a^ 1 {2{z eil v N , x ) NJ -\\v NA \\\ x - IKMl 2 ^} >cr ] 

Unfortunately we are not able to prove this lower bound for general V and W 
but only for unpinned systems (W — 0) in d — 1. We note r T = — g T . Then 
the p's and the r's are independent variables under the Gibbs measure (-)n.t- Some 
conditions on V are imposed 
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Var(V"(r ))<C, Var(V'(r )) < C 

with C > and Var(F(ro)) is the variance (independent of AO of F (ro) under the 
Gibbs measure (-Lt. 

Proposition 2. Under the conditions above there exists a positive constant C' inde- 
pendent ofN and A such that 

C <N- l (3 ei ,(X-J? N )- l 3 ei ) N , T 

Proof. We choose 

v N ,x = -aY,Px(V'(r x ) +V'(r x -i)) = 2a3 ei 

X 

with a > we will precize later. We have 

^ N v N ,x = -a^pl (V"(r x ^)-V"(r x )) 

X 

Let GxAf(z) the solution on T# of 

(A -yA)G N:X = 8o(-) 
Observe that .^{p 2 ) = A(p x ) then we get 

(A - yJ^rVn^) = -aY J G N , l {x-z)p 2 z {y"{r x -y)-V"(r x )) 

x,z 

It follows that 

W N v N , x \l lA = a 2 £ G x , N (x-z)( P y z (V"(r y ^)-V"(r y )) (V"(r x ^)-V"(r x ))) N>T 

x,y,z 

Observe now the Gibbs measure (-)n,t is product. Then an easy computation shows 

1 1 *favn \\- lt x = N ° 2 T Var ( y " fo)) ( 4 G N . X ) (0) 
By discrete Fourier transform we have 

1 _ Asm 2 (nk IN) , 

(4G«)(o) = ttE . - . 2 , L ^y 

A 7 A + 4y sin" (tor /A 7 ) 

On the other hand we have 

2(-Jc[ ,~Vn,x)n,t ~ ||viv,A,|li X 
= 4fl(3 e , ,3 ei )w,r -4a 2 (a e , , (A - y^v)^ ei )iv,r 
= 4a(3 ei )a,t - 4a 2 (A + 7) (3 ei , 3 ei )w,r 
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because S^nZ, 



Zei - By a simple computation we have 



C| IN.T — 



4NTVar(V' '(r )) 



(12) 



If a is sufficiently small we get 



N- d (3 ei ,{x-^ N y l z ei ) N . T >c 



with C > independent of N and A . 



4.3 Disordered harmonic chain 

In this subsection we are interested in the effect of disorder on the thermal con- 
ductivity properties. The simplest way to introduce randomness is to assume that 
masses of atoms can vary from site to site according to a random sequence. We first 
review basic facts for deterministic chains with random masses. As it is well known, 
the presence of disorder generally induces localization of the normal modes and one 
can expect to have a perfect thermal insulators (Kn(T) — > 0). The only analytically 
tractable model is the one dimensional disordered harmonic chain (DHC). Surpris- 
ingly the behavior of the thermal conductivity depends on boundary conditions and 
on the properties of the thermostats (IE), J22)). This curious phenomenon has been 
studied in 1141 (see also f2T| ) in a more general setting and it turns out that "the 
exponent [of Kj\r] depends not only on the properties of the disordered chain itself, 
but also on the spectral properties of the heat baths. For special choices of baths 
one gets the "Fourier behavior" ". If we add a pinning potential in the DHC, Kn(T) 
becomes exponentially small in N. 

Recently, Dhar and Lebowitz ( lfl61P were interested in the effect of both disorder 
and anharmonicity. The conclusions of their numerical simulations are that the in- 
troduction of a small amount of phonon-phonon interactions in the DHC leads to a 
positive finite thermal conductivity. 

We consider now the (a, v) -harmonic chain with random masses and energy 
conserving noise. The noise should simulate in some sense nonlinearity effects and 
in view of the numerical simulations of 0161 one would expect the model to become 
a normal conductor : Kjv(T) — * k(T) with k(T) finite and positive. We are not able to 
obtain interesting informations for k(T) but only for K GK (T). Hence the behavior of 
the thermal conductivity is studied in the linear response theory framework by using 
the Green-Kubo formula. Behavior of the conductivity defined through Green-Kubo 
formula has not been studied for DHC. It would be interesting to know what is the 
order of divergence of the latter. For the perturbed DHC we obtain uniform finite 
positive lower and upper bounds for the d dimensional finite volume Green-Kubo 
formula of the thermal conductivity with or without pinning (Theorem[3]l so that the 
thermal conductivity is always finite and positive. In particular it shows the presence 
of the noise is sufficient to destroy localization of eigen-functions in pinned DHC. 
Linear response approach avoids the difficulty to deal with a nonequilibrium setting 
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where effects of spectral properties of heat baths could add difficulties as in the case 
of purely DHC. In the nonequilibrium setting, we expect that since the Green-Kubo 
formula for the thermal conductivity of the perturbed DHC remains finite, it will not 
depend on the boundaries. 

Theorem 3. There exists a positive constant C > independent of X andN such 
that 

CT 1 < liminfliminf f e- kt N~ d E NT [Z e , (t), Ze, (0)]dt 

A^O AT-m» Jo 1 

< limsuplimsup f e~ Xt N~ d E NiT [Z ei (f), Z ei (0)]dt < C 

X^O N^oo JO 

Proof. The proof is similar to what is presented in subsection l4.2l for the anharmonic 
case. 

A priori the Green-Kubo formula K GK {T) depends on the particular realization 
of the masses {m x }. But a formal argument (J6|) suggests that if masses are dis- 
tributed according to some stationary ergodic probability measure P* then K GK {T) 
depends only on the statistics of the masses and not on the particular realization of 
the disorder. One can write formally an infinite Green-Kubo homogenized formula 
Khom. (T) obtained by averaging over the masses. In order to define JQj 0m . (T) one has 
to consider the infinite volume dynamics with generator _Sf given by 

I {dp,*? ■ d qx - d qx W ■ d Px ) + 7 E E (X l x i +ek ) 2 
xez d ij-k=ixez d 

and 

P 2 

JT= E —+{aA-v)q x -q x 

Observe that the sums are taken over 7L d . One can show that the dynamics with 
generator Jz? is well defined. Then we have 



A — >U 



o 



dte-^^E T [jl ei (t)^jl ei 



where T- denotes the shift on the configuration space and E r is the expectation cor- 
responding to the infinite dynamics starting from the infinite volume Gibbs measure 
with temperature T. 

Theorem 4. Assume that {m x } is stationary under P* and there are positive 
constants m and m such that 

P*(m<m x <m) = 1 

The Hamiltonian contribution to the homogenized Green-Kubo formula for the ther- 
mal conductivity »c££ (T) — y/d 
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^(rj-y/^HmE 




zez d 



exists, is positive and finite. 

Proof. The proof is based on functional analysis arguments (||4), IfTO I, [6 1). 



5 Energy-momentum conserving model 

We investigate the same problems as in the previous section but for the chain per- 
turbed by the energy-momentum conserving noise. Momentum is then conserved by 
the total dynamics if it is conserved by the Hamiltonian dynamics which is the case 
if and only if the system is unpinned. Dramatic consequences in low dimensional 
systems appear. In fact we prove that for unpinned systems with harmonic interac- 
tions, thermal conductivity is infinite in 1 and 2 dimensions, while is finite for d > 3 
or for pinned systems. So thermal conductivity in our model behaves qualitatively 
like in a deterministic nonlinear system, i.e. these stochastic interactions reproduce 
some of the features of the nonlinear deterministic hamiltonian interactions. 

For anharmonic systems, even with the stochastic noise, we are not able to prove 
the existence of thermal conductivity (finite or infinite). If the dimension d is greater 
than 3 and the system is pinned, we get a uniform bound on the finite size system 
conductivity. For low dimensional pinned systems (d = 1 , 2), we can show the con- 
ductivity is finite if the interaction potential is quadratic and the pinning is generic. 
For the unpinned system we have to assume that the interaction between nearest- 
neighbor particles is strictly convex and quadratically bounded at infinity. This be- 
cause we need some informations on the spatial decay of correlations in the sta- 
tionary equilibrium measure, that decay slowly in unpinned system. In this case, we 
prove the conductivity is finite in dimension d > 3 and we obtain upper bounds in 
the size N of the system of the form \/N in d = 1 and (\ogN) 2 in d = 2 (see Theorem 



5.1 Homogenous harmonic chain 

Consider first the homogenous (a, v)-harmonic case (]]) with all masses equal to 1. 
In we obtained the following theorems. 

Theorem 5. ^2]/ In the (a, v)-harmonic case (0, the limit defining K GK (T) exists. 
It is finite if d > 3 or if the on-site harmonic potential is present (v > 0), and is 
infinite in the other cases. When finite K GK (T) is independent of T and the following 
formula holds 
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where 0){k) is the dispertion relation 



) 



1/2 



a>(k) = v + 4a sin 2 (nk j ) 



and 



W(k) 



( 



8lf = isin 2 (^'), if d>2 
4/3 sin 2 (7TA:)(l+2cos 2 (7rA:)), 



if d=l 



Consequently in the unpinned harmonic cases in dimension d = 1 and 2, the con- 
ductivity of our model diverges. We have the following behavior for the truncated 
Green-Kubo formula. 

Theorem 6. ft2§ In the harmonic case, ifW = 0: 



In all other cases KjJ is bounded in N and converges to K GK {T). 

The interest of these theorems is that they show energy-momentum conserving 
model reproduces the expected behavior of purely Hamiltonian nonlinear chains, 
meaning an infinite conductivity for low dimensional unpinned systems and a finite 
conductivity otherwise. Several microscopic stochastic models have been proposed 
in the past but none of them has this property with respect to momentum conserva- 
tion and dimension. 

The strategy of the proof is very similar to the proof given for the energy conserv- 
ing model. Recall (0. Then the problem is reduced to solve the resolvent equation 



Solving explicitly such an equation is in general very difficult. The key property of 
_S^V is that for quadratic potentials V and W it sends a polynomial function of p's 
and q's, of degree A: to a polynomial function of the same degree. Since Ze\ is of 
degree 2, u N i has to be searched in the smaller space of polynomial functions of 
degree 2. Moreover the translation invariance of 3 Cl implies that u N i is in the form 

U N.X = fN,X i x - y) Px ■ Py + Y, h N,xi x -y) <7-v ' <ly + £ 8N, X (* ~ ?) Px ' <7v 



where f N XiiN,x an d n N.i are functions from Tf, to R. Then we obtain linear equa- 
tions for fn,x,gN,x and /i^A. It turns out that f N x =fyvA = and that g N i is simply 
related to the Green function of a symmetric simple random walk on T^. Then one 
computes easily 



1. K GK {T)^N X I 2 ifd= 1, 

2. K$ k \t) ~ logA? ifd = 2. 




and we get the result. 
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For the anharmonic chain we are not able to solve the resolvent equation 

%<u N: x — J^nUn a = 3pi 

Nevertheless we can use the same strategy as the one explained in subsection 14.21 
We introduce a variational formula to estimate by above 

( u N,l i -Jei )n,T 

Extra assumptions on the potentials V and W assuring a uniform control on the 
canonical static correlations have to be done. We have this control as soon as V is 
strictly convex. In the pinned case W > 0, this control is "morally" valid as soon 
as the infinite volume Gibbs measure is unique. Exact assumptions are given in 
|9l , theorem 3.1 and theorem 3.2. Hence "general anharmonic case" will refer to 
potentials V and W such that this control is valid. 

Theorem 7. ^2]/ Consider the general anharmonic case. There exists a constant C 
(depending on the temperature T) such that 

• For d > 3, 

1. either Wj > is general 

2. or ifWj = and < c_ < V" <C+< °°for any j, 
then 

*N K (T) < C. 

• For d = 2, ifWj = and < c_ < V'j < C+ < °°for any j, then 

*N K {T) < C(logA0 2 . 

• For d = 1, ifWj = and < c_ < V'J < C+ < «, then 

K$ K (T) < CVN. 

• Moreover, in any dimension, if Vj are quadratic and Wj > is general then 
4 K (T)<C. 



6 Open questions 

In this section we mention several open questions. From the (probably) easier to the 
most difficult we have: 
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• Prove © in the li-dimensional linear pinned or unpinned case. It should be also 
valid in the nonlinear case. 

• Consider the energy-momentum conserving model with linear interactions. We 
have studied the thermal conductivity in the Green-Kubo framework. Consider 
the same system but in contact with thermal baths at different temperatures. Can 
you say something about Kfr(T) ? We expect it diverges like \f~N (resp. logAO 
for d = 1 (resp. d = 2) in the unpinned case and converges in the other cases 
to a finite positive constant. In fact a similar unpinned one-dimensional model 
has been considered in lfT3l . |[T8l and numerical and analytical (but not rigorous) 
results confirm this conjecture. 

• Consider the (a, v) -harmonic case with random masses and energy conserving 
noise. Can you say something about the temperature dependance of ?cf™\ ? Can 
you prove the almost sure convergence (w.r.t. disorder) of the Green-Kubo for- 
mula to the homogenized Green-Kubo formula fc££ ? 

• Consider the (a, v)-harmonic case with random masses and energy-momentum 
conserving noise. We have still upper bounds for the truncated Green-Kubo for- 
mula similar to the one obtained for the homogenous nonlinear chain in theorem 
|7] Can you obtain lower bounds? In particular do you have a positive conductiv- 
ity for pinned chains ? 

• Can you prove convergence of the Green-Kubo formula for nonlinear chains with 
energy conserving noise ? 

• Consider an homogenous nonlinear chain with energy-momentum conserving 
noise. Can you obtain (non trivial) lower bounds for the truncated Green-Kubo 
formula ? Can you prove convergence of the Green-Kubo formula in dimension 
d > 3 or for pinned systems ? Interesting and surprising numerical simulations 
are provided in O . 

• Consider an homogenous nonlinear chains with energy-momentum conserving 
noise in contact with two thermal baths at different temperatures. Can you say 
something about Kjv(r) ? 

• Consider a nonlinear chain with energy-momentum conserving noise and random 
masses. Can you say anything about k(T) or K GK (T) ? 
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